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FOREWORD 


This  report  covers  in-house  research  conducted  by  the  Air  Force  Flight 
Dynamics  Laboratory,  Vehicle  Equipment  Division,  Environmental  Control 
Branch.  This  study  effort  is  in  support  of  ProjectNo.  1309,  "Environmental 
Interactions, "  Task  No.  130904,  "Vibration  Technology, "  and  was  conducted 
under  the  direction  of  Robert  W.  Sevy  and  David  A.  Earls.  Experimental  as¬ 
sistance  was  provided  by  Messrs.  E.  Ruddell  and  F.  Stahl;  data  analysis  and 
reduction  by  Messrs.  E.  Schell  and  H.  Kinum.  This  report  represents  the 
work  done  during  the  period  June  1965  through  March  1968. 

This  technical  report  has  been  reviewed  and  is  approved. 


WILLIAM  C.  SAVAGE  ^ 

Chief,  Environmental  Control  Branch 
Vehicle  Equipment  Division 
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ABSTRACT 

This  study  is  primarily  concerned  with  the  prediction  of  internal  vibration 
levels  in  an  equipment  package  using  statistical  energy  methods.  Modal  densi¬ 
ties,  loss  and  coupling  factors  were  determined  theoretically  and  experimen¬ 
tally.  Predicted  vibration  levels  were  derived  and  compared  to  experimental 
results.  The  energy  coupling  factor  relating  the  box  cover  to  the  internal 
component  responses  was  derived  and  compared  favorably  to  experimental 
values.  A  generalized  coupling  factor  is  included  for  future  applications.  It 
is  demonstrated  that  statistical  energy  methods  offer  a  new  and  powerful 
method  by  which  the  average  internal  vibration  levels  of  equipments  may  be 
predicted  when  the  vibration  field  of  the  equipment  cover  is  known  or  can  be 
estimated. 

This  abstract  is  subject  to  special  export  controls  and  each  transmittal  to 
foreign  governments  or  foreign  nationals  may  be  made  only  with  prior  approval 
of  the  Air  Force  Flight  Dynamics  Laboratory  (FDFE),  Wright -Patterson  Air 
Force  Base,  Ohio  45433. 
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SYMBOLS  (CONTD) 

=  mass  density  of  air 
=  velocity  of  sound  in  air  (in/ sec) 

=  frequency  (Hz) 

-  circular  frequency  (rads/sec) 

=  bandwidth  (rads/sec) 

2 

=  area  of  box  or  boards  (in  ) 
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1 

V 

=  modal  resonance  frequency  (Hz) 

APPENDIXES  I,  II,  AND  III 

=  spring  constant  of  stripped  circuit  boards  (lbs/in) 

=  spring  constant  of  fully  loaded  and  mounted  boards  (lbs/in) 

=  modulus  of  elasticity  (stripped  board) 

in 

=  modulus  of  elasticity  (boards  loaded  and  mounted) 

=  weight  of  board  without  components 
=  weight  of  board  with  components  and  stiffeners 
=  static  force  at  ^  (lbs) 

=  midpoint  of  circuit  boards  (in) 
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A  =  deflection  of  board  at  midpoint  (in) 

b  =  board  width  (in) 
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SECTION  I 

«  INTRODUCTION 

Intense  sound  and  vibration  fields  are  generated  by  the  propulsion  systems 
of  aircraft,  missiles,  and  by  high  speed  flight  regimes  resulting  in  air  pressure 
fluctuations  which  in  turn,  couple  to  structural  surfaces.  This  energy  couples  to 
and  through  equipments  to  produce  complex  responses  which  may  result  in  mal¬ 
function  of  internal  systems  and/or  in  fatigue  failures.  In  order  to  construct  an 
item  of  equipment  to  withstand  the  acoustic  and  vibration  environment,  a  pre¬ 
diction  process  is  required  that  permits  the  designer  to  estimate  the  response 
magnitudes  of  the  internal  levels  over  a  broad  frequency  range. 

The  application  of  statistical  energy  methods  to  flight  vehicle  equipment 
provides  a  promising  approach  for  predicting  vibration  levels  at  various  regions 
of  the  equipment  structure  and  substructure.  By  use  of  power  flow  analogies, 
the  coupling  behavior  between  the  vibration  source  and  receiver  may  be  de¬ 
scribed  and  then  utilized  to  define  the  time  and  spatially  averaged  vibration 
magnitudes  which  result  from  an  external,  coupled  excitation  representing  an 
acoustic  field  or  another  vibration  source. 

This  study  effort  consists  of  six  sections.  The  first  section  contains  a 
short  introduction  and  the  second  section  gives  the  approach  to  the  study  effort. 
The  third  section  is  a  study  of  the  modal  densities  of  the  box  cover  and  the 
boards  and  is  concerned  with  the  development  of  a  general  equation  by  which 
the  mode  densities  of  printed  circuit  boards  may  be  adjusted  for  mass  loading 
effects  and  stiffness  influences  on  the  modal  count  due  to  the  presence  of 
mounted  components  and  due  to  the  addition  of  support  stringers.  The  fourth 
section  is  concerned  with  the  measurement,  and  derivation  of  the  applicable 

2 

loss  factors  that  are  required  in  order  to  obtain  the  energy  coupling  factor, /x^  . 

The  fifth  section  proposes  the  prediction  of  the  vibration  response  of  the  vibration 

response  of  the  boards  resulting  from  a  known  vibration  field  of  the  box  cover. 

The  predicted  results  are  compared  with  a  set  of  laboratory  measurements. 

2 

The  experimental  /x  1  is  also  obtained  and  compared  to  the  derived  form.  A 
2  1 

smoothed  /x^  is  synthesized  for  tentative,  future  applications.  Finally,  the 
sixth  section  consists  of  conclusions  and  recommendations  for  future  studies. 
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SECTION  II 
APPROACH 


Consider  a  two  element  energy  sharing  structure  which  might  represent, 
one,  the  outer  cover  of  an  equipment  box  (symbolized  by  subscript  1),  and, 
two,  the  inner  connecting  elements  consisting  of  printed  circuit  boards  (sym¬ 
bolized  by  subscript  2). 

2 

From  Reference  1  we  obtain  the  energy  relation  TT2  =  tt1’  where 

T_0  and  Tr_1  are  the  total  kinetic  energies  of  corresponding  mode  sets  in 
1  ^  li  2 

structures  1  and  2  and  where  represents  the  energy  coupling  factor.  The 

factor,  ,  is  equal  to  V2i^  \\  +  ’’V  where  \\  is  the  couPlinS  loss 

factor  arising  from  the  mechanical  joints  between  structures  1  and  2  and  where 

17  is  the  true  loss  factor  of  structure  2  taken  alone  (removed  from  structure  1). 
z 

m1  —  2  m2  _  2 

Again,  in  Reference  1  we  note  that  T1  =  >  andT^^j^-  <V2  > 

1  2 

where  N1  =  nx  (cu)  Acu  and  Ng  =  n2  (cu)  Acu.  N ±  is  the  total  number  of  modes  in 

energy  set  one  and  N2  is  the  total  number  of  modes  in  energy  set  two  and  Acu 

is  the  set  bandwidth.  T0  and  T..  are  the  average,  single  modal  energies.  Sub- 

2  ^  ^ 

stituting  in  T2  =  fj.  ^  ^  to  solve  for  the  mean  squared  velocity  of  the  modal 
response,  we  have: 

<  V2  >  _  n2  (cu)  m,  2 

- — - -  •  /x,  U) 

<  V|  >  Dj  (cu)  m2 

Where  n(<u)  is  the  structural  mode  density,  m^  represents  the  structural 
mass  and  <V‘‘J>  is  the  space  and  time  averaged  mean  square  velocity.  For  the 
case  of  random  vibration  taken  over  a  mode  ensemble  of  bandwidth  Acu  ,  we  set 


<  V2>cu2 
Acu 


<S  (<u)> 


which  is  the  spatially  averaged  acceleration  power  spectral  density.  Equation  1 
now  expresses  the  response  of  the  printed  circuit  boards  as  follows: 


<S(cu)2>  =  <  S  (cu) ,  > 


n2(cu)m,  ; 
n,  (cu)m2 


(2) 
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where: 


2 

Ml  = 


"721 

V2\  1?2 


(3) 


Similarly,  by  equating  the  power  flow  between  an  acoustic  volume  enclosing  an 
equipment  box  one  may  describe  the  resultant  vibration  field  of  the  box  cover. 
From  Reference  1  we  have, 


<  P2> 


CQ  ni(aj) 

PQ  m' 


(4) 


o 

Where  <  a  >  is  the  same  spatial  and  time  averaged  mean  squared  acceleration 

of  the  box  cover  resulting  from  a  time  and  spatially  averaged  mean  squared 

—2 

acoustic  pressure  <p  >  .  The  terms  c  and  p  are  the  velocity  of  sound  and 

a  a  2 

the  density  of  air.  The  term  p  is  more  complex  than  p^  reflecting  the  coupling 
from  box  cover  to  board  as  well  as  the  energy  coupling  of  the  reverberant  room 
to  the  box.  It  is  not  intended  to  consider  Equation  4  further  in  this  study  other 
than  to  show  that  by  the  use  of  energy  equations,  we  may  be  able  to  predict 
vibration  levels  of  the  internal  elements  of  an  equipment  package  from  a  directly 
exdited  acoustic  volume.  In  this  study  we  are  chiefly  concerned  with  the  pre¬ 
diction  of  the  vibration  level  of  the  equipment  boards  when  the  vibration  level 
of  the  box  surface  is  known  or  can  be  estimated. 


Returning  to  Equation  2  and  the  equipment  box,  we  see  that  the  prediction 

of  <S(u>)  >  depends  upon  the  determination  of  the  modal  densities  of  both  box 
2  2 
cover  and  boards,  their  masses,  and  the  energy  coupling  factor,  p^  .  Since 

structural  masses  are  usually  known,  or  can  be  estimated,  it  remains  to  de¬ 
termine  n^(cu),  n^t^uj.and  p,  if  the  prediction  process  is  to  be  utilized.  The 
details  of  this  study  involves  the  measurement  and  derivation  of  the  modal 

densities,  together  with  a  description  of  a  variety  of  loss  factors  whose  inter- 

2 

relationships  allow  us  to  derive  p^ .  These  properties  are  then  applied  to  an 
aircraft  equipment  which,  reduced  to  its  simplest  model,  is  presumed  to  con¬ 
sist  of  an  array  of  printed  circuit  boards  enclosed  by  a  sheet  metal  cover. 

The  equipment  selected  is  an  electronic  equipment  package  and  is  part  of 
the  gunfire  alarm  and  gun-fire  detection  system  for  fighter  aircraft.  It  contains 
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transistorized  circuitry  with  ten  plug-in  printed  circuit  boards  stacked  in  two 
vertical  columns.  The  box  size  is  15  3/8"  x  9  7/8"  x  7  5/8"  and  weighs  18 
pounds.  Figure  1  shows  the  item  with  the  cover  removed. 
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Figure  1.  Equipment  Box  With  Side  Cover  Removed 
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SECTION  III 
MODAL  DENSITY 


1.  DEFINITION 

The  theoretical  expression  (Equation  3)  for  the  modal  density  of  a  simply 
supported  isotropic  thin  plate  is: 


n  (cu) 


AP  _ P 

47rKCL  2  y<uKCL 


(5) 


where,  Ap  is  the  surface  area  (taken  on  one  side),  P  is  the  surface  perimeter, 
CL  is  the  longitudinal  wave  velocity,  t  is  the  plate  thickness;  and  K  is  the  radius 
of  gyration  and  is  equal  to 


It  has  been  shown  in  past  studies  that  the  second  term  is  negligible  (over  most 
of  the  frequency  range)  and  may  be  dropped.  Thus,  the  mode  density  becomes: 


n  (cu) 


n  (f) 


47rKC[_ 

APyr 


(6) 

(7) 


The  next  step  is  to  examine  the  tractability  of  the  modal  density  expression 
when  applied  to  loaded  printed  circuit  boards  and  actual  box  covers. 


2.  MODAL  MEASUREMENTS  OF  BOARDS  A  AND  B 


Two  printed  circuit  boards  were  removed  from  the  box  and  measurements 
were  made  of  their  modal  densities.  One  plug  in  board  (Figure  2)  approxi¬ 
mately  8.5"  x  1/16"  was  covered  by  resistors,  capacitors,  and  transistors. 
The  second  board  is  approximately  8.  5"  x  6"  x  1/16"  and  is  heavily  loaded 
with  shift  register  cores  along  with  resistor  elements.  For  ease  of  identifica¬ 
tion,  we  shall  hereafter  refer  to  each  board  type  as  A  and  B  respectively. 
Both  boards  are  plug-ins  and  are  secured  at  the  end  opposite  the  tongue  with 
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Figure  2.  Plug-In  Boards,  Types  A  and  B 
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two  screws  located  in  5/16"  square  hollow  aluminum  rods.  Boards  A  and  B 
are  edge  bound  by  1/16"  aluminum  angle  strips  that  are  attached  to  the  board 
edges  with  widely  spaced  rivets. 

A  jig  was  constructed  to  hold  the  boards,  to  facilitate  modal  density 
measurements  of  the  boards  when  removed  from  the  box,  and  to  allow  close 
inspection  of  the  modal  patterns  (Figure  3).  Mounting  holes  were  located  and 
a  mating  plug-in  receptacle  was  positioned  to  simulate  the  in  situ  mounting 
conditions  with  the  important  exception  that  the  jig  structure  is  sufficiently 
rigid  as  to  be  considered  relatively  lossless  beyond  the  boundaries  of  the  two 
mounting  screws  and  the  plug-in  receptacle.  Accelerometers  (Endevco  2226) 
were  mounted  with  Eastman  910  cement.  Several  locations  per  board  were 
recorded.  The  amplified  signal  (Endevco  2614B)  was  coupled  into  a  Spectral 
Dynamics  SDS-101A  analyzer  using  a  2  CPS  crystal  bandwidth  chosen  for  good 
resolution  of  the  responses.  The  filtered  output  was  coupled  to  a  Hewlett- 
Packard  logarithmic  converter  (Model  60D);  thence,  to  the  Y  axis  of  11"  x  17" 
Electro  instruments  Model  500  X-Y  recorder.  The  X  or  frequency  axis  was 
supplied  from  a  Hewlett-Packard  500c  frequency  converter  to  allow  expansion 
of  the  frequency  axis  and  thus  augmenting  detailed  resolution  of  modal  re¬ 
sponses  and  also  providing  a  more  accurate  measurement  of  the  loss  factor 
when  using  the  half  power  method. 

A  Ling  Electronics  A-300B  shaker  was  used.  The  board  was  subjected  to 
lg  peak  sinusoidal  acceleration;  slowly  swept  from  20  to  2000  CPS.  The  board 
was  checked  over  the  frequency  range  for  linearity  of  response  by  increasing 
the  input  level  and  noting  a  proportionate  increase  in  the  accelerometer 
response. 

The  responses  were  recorded  on  11"  x  17"  graph  paper.  Three  records 
per  total  sweep  were  recorded:  25  to  100  CPS,  100  to  1000  CPS,  and  1000  to 
2000  CPS. 

3.  MODAL  DENSITY  OF  BOARDS  A  AND  B 

The  accumulative  mode  count  for  boards  A  and  B  are  shown  in  Figures 
4  and  5.  A  straight  line  was  drawn  through  the  points.  The  slope  of  the  curve 
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Figure  3 


Test  Jig  Used  for  External  Modal  Measurements 
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Figure  4.  Cumulative  Mode  Count  for  Printed  Circuit  Board  A,  Externally 
Mounted 
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Figure  5.  Cumulative  Mode  Count  for  Printed  Circuit  Board  B,  Externally 
Mounted 
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represents  the  observed  modal  density  n(f)  .  The  dashed  line  represents  the 
theoretical  n(f)  (of  the  stripped  boards)  using  Equation  7. 

It  is  seen  that  the  theoretical  values  exceed  considerably  the  measured 
value  s  -  approximately  50  percent  for  board  A  and  20  percent  for  board  B. 


4.  CORRECTIVE  STEPS 


Assuming  that  the  difference  between  measured  and  calculated  values  is  due 
primarily  to  the  mass  loading  of  the  mounted  components  and  the  increased 
stiffness  due  to  components,  stiffeners,  and  other  mounting  stiffness  contribu¬ 
tions,  we  proceed  to  attempt  a  correction. 


Recalling  that. 


n  (f) 


AP  V^3~ 
CL  t 


Where, 


n  (f)  = 

AP  = 
t  = 


mode  density  (modes/CPS) 
the  boord  area  (in.^) 
board  thickness  (in.) 


yiE 

p 


the  longitudinal  wave  velocity  (in. /sec) 


We  note  that  C^  describes  the  stiffness  and  mass  properties.  A  new  board  is 
imagined  having  a  distributed  mass  and  stiffness  equivalent  to  our  fully  loaded 
boards  mounted  in  situ  (the  jig  mounting  stiffness  is  considered  to  approximate 
mounted  conditions  in  the  box).  For  the  stripped  boards  we  let: 


cL  t 

and  allow:: 


n(f') 


ap  yr 

CL  * 
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The  egression,  n(f  )',  represents  the  mode  density  of  the  fully  loaded  boards 
having  an  adjusted  longitudinal  wave  velocity  reflecting  the  stiffness  and 
mass  contributions.  For  the  derivation  of  this  expression  see  Appendixes  I 
and  II.  Dividing  one  equation  by  the  other  and  solving  for  n(f  ),  we  obtain, 


where  c'  is  termed  the  adjusted,  or  pseudo  velocity. 

L 

For  board  A,  (computations  in  Appendixes  I  and  II) 

n(f)  is  approximately  18.  2  modes/1000  CPS 
and  n(f  )'  is  12. 1  modes/1000  CPS 

For  board  B ,  (Appendixes  I  and  II) 

n(f)  is  approximately  20  modes/1000  CPS 
and  n(f  )’  =  17.  7  modes/1000  CPS 

The  corrected  mode  densities  n(f  )'  are  compared  to  the  observed  values  n(f)Q 
and  the  computed  stripped  board  values  n(f).  The  plots  of  n(f)'  is  also  included  in 
Figures  4  and  5. 

BOARD  A  BOARD 

n(f)  =  18.  2  modes/1000  CPS  n(f)  =  20.  0  modes/1000  CPS 

n(f  )'  =  12. 1  modes/1000  CPS  n(f  )'  =  17.  7  modes/1000  CPS 

n(f)Q  =  12.  0  modes/1000  CPS  n(f)Q  =  17.  3  modes/1000  CPS 

a.  Results 

In  both  cases  the  corrected  value  of  the  modal  densities  are  in  much 
closer  agreement  with  the  observed  estimates. 

It  does  appear  that  acceptable  engineering  estimates  of  the  modal  density 
of  heavily  laden  panels  (including  stiffeners)  can  be  realized  if  the  longitudinal 
wave  velocity  of  the  panel  (taken  alone)  can  be  adjusted  by  the  estimated  or 
measured  stiffness  and  mass  factors  of  the  actual  board. 
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5.  MODAL  OVERLAP  ANALYSIS 


We  seek  reasonable  assurances  that  n(f)Q,  the  observed  modal  density,  is 
approximately  equal  to  the  true  modal  density,  n(^rue*  That  is,  within  the 
frequenc:y  range  of  this  study  was  there  sufficient  modal  overlap  to  conceal 
more  modes,  than  were  observed? 


If  statistically,  the  average  bandwidths  of  modes  (e.  g.  bandwidth  at  the 
3  db  down  points)  approaches  the  modal  separation,  ,  then  modal  overlap 
should  occur  and  n(f)Q  should  be  less  than  n(f)j.rue*  Restating,  when  =  Ab 
we  should  expect  that  on  the  average  n(f)Q  would  be  approximately  one-half  of 
n(f)true  —  1x1  any  case,  less  than  the  ^true*  This  is  illustrated  in  Figure  6. 

We  examine  the  region  of  the  frequency  spectrum  at  which  this  might  occur. 
From  Reference  2 


let  the  modal  separation,  Af,  be  equal  to  the  modal  bandwidth,  Ab 


where, 


Af 


I 

n  ^  ^True 


Ab 


At>  =  -£-ond  Q  =  -i- 

substituting  Ab  into  Equation  8  and  for  Q  we  obtain; 

f  =  _ ! _ 

n  17  ,  n  (f ). 

'  ’  true 

where  f  is  the  frequency  above  which  modal  blending  is  pronoimced. 
n 


(8) 


(9) 


a.  Results 


A  review  of  the  Ab's  obtained  from  the  vibration  records  show  that  over 
the  vibration  range  the  modal  bandwidths  (as  measured  by  the  analyzer)  were 
on  the  whole  less  than  the  observed  or  computed  modal  separation. 
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Figure  6.  Modal  Separation  and  Blending 
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For  example: 


Recalling  that,  AfQ 
bandwidths 


n(% 


;  we  compare 


modal  separation  to  response 


Af  Ab 

o  o 


Board  A 

83  CPS 

from  6  to  60  CPS 

Board  B 

59  CPS 

from  4  to  58  CPS 

These  results  indicate  that  relatively  little  blending  of  modes  occur  -  and  that 
n(f)  represents  a  reasonable  estimate.  The  records  tend  to  confirm  this.  On 
the  whole,  the  modes  are  distinct  and  identifiable  up  to  2000  CPS.  To  estimate 
the  frequency  above  which  we  would  expect  substantial  overlap,  we  return  to 
Equation  9 


I 

17  ,  n  (f )  true 


set  tinqn(f)  =n(f). 

true  0 


For  board  A 


and  board  B 


n(f)Q  =  0.012  modes/CPS 


n(f)0  =  0.017  modes/CPS 

It  will  be  shown  in  Section  IV  and  in  Figure  13,  that  the  loss  factors  of  both 
boards  at  the  higher  frequencies  (2000  CPS)  tend  to  be  independent  of  the  fre- 

_2 

quency  and  approximately  asymptotic.  In  our  case  it  is  approximately  2.  0  x  10  . 


Then, 


And, 


f.  =  - — - I03  ~  4200  CPS 

nA  ( 2.0 )  ( 12  ) 


100 

(2.0)117) 


10' 


~  2900  CPS 


Of  course,  some  modal  overlap  occurs  even  below  2000  CPS.  Further,  the 
overlap  criterion  of  Equation  9  requires  that  the  response  magnitudes  be  suf¬ 
ficiently  high  to  obtrude  from  the  skirts  of  adjacent  responses,  thus  modal 
overlap  is  also  related  to  response  magnitude. 
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b.  Conclusions 

Observed  modal  counts  should,  therefore,  be  considered  as  minimum 
values.  It  may  be  significant  that  in  both  cases  the  adjusted  values  of  the  modal 
density  of  the  boards  remain  somewhat  greater  than  the  observed  counts  and 
that  this  difference  represents  occasional  modal  blending. 

6.  MODAL  DENSITY  OF  OUTER  PANELS 

Modal  responses  of  the  panels  were  recorded  through  a  2  CPS  filter.  Three 
different  locations  on  each  face  of  the  box  were  recorded.  These  mode  counts 
were  averaged  per  panel  and  then  averaged  for  six  panels  to  provide  a  total 
mean  mode  count  for  the  box  cover.  In  these  experiments  one  small  high  fre¬ 
quency  shaker  (MB  1250)  became  available.  The  shaker  was  attached  to  a  side 
panel  of  the  box  and  remained  in  this  location  throughout  the  recording. 

a.  Results 

Cumulative  modes  from  the  small  shaker  experiments  were  plotted  in 
Figures  7  and  8.  For  graphical  clarity  the  data  points  of  Figure  7  are  omitted. 
The  cumulative  mode  plot  is  shown  as  a  smoothed  curve  with  a  positively  de¬ 
creasing  slope  above,  approximately  300  CPS.  Figure  8  represents  the  expansion 
of  Figure  7  in  the  low  frequency  range.  The  details  of  the  expansion  emphasize 
the  difficulty  of  determining  the  modal  density  when  modal  blending  is  encoun¬ 
tered.  Below  300  CPS,  the  cumulative  mode  history  is  reasonably  linear  and 
behaves  acceptably.  But  above  this  range,  the  observed  slope  n(f)o  decreases; 
giving  slope  values  approximately  one-half  that  of  the  initial  observation. 

Examination  of  the  vibration  records  showed  a  distinct  trend  of  modal 
overlap  above  approximately  500  CPS  —  this  is  illustrated  by  samples  (a),  (b), 
and  (c)  of  Figure  9.  The  certainty  of  questionable  modes,  indicated  with  arrows 
on  (b)  and  (c),  was  verified  by  inspection  of  other  accelerometer  records  which 
usually  presented  these  obscured  modes  in  much  more  distinct  form.  In  fact, 
the  multiple  record  method  was  used  extensively  when  such  doubts  arose. 
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Figure  7.  Mode  Density  of  Cover  Panels  Showing  Effect  of  Modal  Overlap 
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Figure  8.  Mode  Density  of  Panels  Below  Overlap  Region 
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Figure  9.  Typical  Modal  Responses  Showing  Increasing  Modal  Overlap  Vs. 
Frequency 
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It  is  seen  then  that  the  decreasing  slope  above  500  CPS  is  due  to  the  modal 
sharing  process  and  that  the  initial  slope  rather  than  the  final  slope  represents 
the  most  accurate  estimate  of  the  modal  density. 


As  in  the  case  of  the  printed  circuit  boards,  one  might  estimate  the  modal 
density  box  cover  from  the  expression, 


n  l  f ) 


cLt 


This  would  be  a  minimum  value  due  to  the  contributions  from  mechanical  cou¬ 
pling  through  and  to  the  printed  circuit  boards.  Thus,  the  actual  mode  density 
should  yield  an  observed  n(f)Q  greater  than  the  theoretical  estimate,  n(f).  We 
investigate  this  returning  to  the  modal  density  expression,  where,  for  aluminum, 

C,  =  200,  000  in/sec,  t  =  panel  thickness  (.  065  inches)  and  A  is  the  area  of  the 
L  P 


panels.  We  have  in  tabulated  form: 

Length  x  Width  (in)  Area  (in2)  n(f ) 

Top  and  bottom  panel  15.  375  x  9.  75  x  2  149.  91  x  2  39.  90 

Side  panel  (1)  15.  375  x  7.  50  115.  32  15.  35 

End  panel  (2)  7.  50  x  9.  75  x  2  73.  13  x  2  19. 46 

Cover  panel  (1)  15.  25  x  7.  50  114.  38  15.  22 

Total  Theoretical  n(f)  =  89.  93  x  10^ 


The  observed  modal  density  of  the  panels  should  then  be  at  least  90  x  10  . 

_3 

Examination  of  the  slope  of  Figure  8  shows  an  n(f)Q  of  105  x  10  in  rather 
close  agreement  with  n(f). 


b.  Conclusions 

The  initial  slope  of  the  accumulative  mode  curve  for  the  box  cover  repre¬ 
sents  the  most  accurate  estimate  of  the  mode  density. 

The  mode  density  of  the  cover  panels  is  in  good  agreement  with  the  theo¬ 
retical  expectation;  being  slightly  less  but  approximately  equal  to  the  observed 
value. 
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7.  SUMMED  MODAL  DENSITY  OF  BOARDS 

Ultimately,  we  are  concerned  with  the  boards  not  as  individual,  discon¬ 
nected  energy  sinks  but  rather  as  a  two-column  parallel  array  of  boards  having 
a  total  mode  density  equal  to  the  sum  of  the  individual  boards  —  this  is  the 
meaning  of  element  two  of  our  chosen  model;  the  two  element  energy  sharing 
structure.  Therefore,  to  obtain  )  of  Equation  2  we  add  the  board  modal 
densities;  observing  that  for  like  boards  we  may  obtain  the  totals  in  proportion 
to  their  areas. 

Summing  the  mode  densities  of  the  boards,  we  have  4  boards  of  type  A  and 
3  of  the  type  B ,  also  three  smaller  boards  similar  to  type  A  and  totaling  1. 45 
times  the  area  of  a  type  A  board. 

n(f)Q  of  A  =  4  x  12.  0  =  48.  0  x  l(f 3 

n(f)Q  of  small  boards  =  1.  45(12.  0)  =  16. 4  x  10-3 

n(f)Q  of  B  =  3  x  17.  3  =  51.  9  x  10~3 

Total  n(f)  =  116.  3  x  10~3 

o 

Here,  we  have  used  the  observed  board  densities  but  we  could  just  as 

easily  have  substituted  the  adjusted  values,  n(f)'.  This  would  result  in  a  total  of 
-3 

119  x  10  —  an  insignificant  difference  representing  about  two  percent. 

In  retrospect,  we  note  that  the  predicted  modal  density  of  the  box  cover 

_3 

is  about  90  x  10  ,  being  approximately  equal  to  that  of  the  board  assembly. 

Later  on  (Sections  IV  and  V)  we  will  set  n^(f)  =  ng(f),  obtaining  small  error 
and  gaining  some  manipulative  advantage  as  a  result.  It  should  be  stressed 
that  it  is  not  necessary  that  this  near  equality  exist  —  it  simply  happens  to  be 
so,  in  this  case. 

8.  MODAL  PATTERNS  OF  BOARDS  A  AND  B 

After  mode  density  measurements  had  been  taken  in  the  test  jig,  each 
board  was  sprinkled  lightly  with  fine  Ottowa  sand  and  modal  patterns  and  their 
frequencies  were  recorded  (Figures  10  and  11). 


22 


AFFDL-TR-69-54 


23 


Figure  10.  Modal  Patterns  of  Printed  Circuit  Board  A 
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Figure  11.  Modal  Patterns  of  Printed  Circuit  Board  B 
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a.  Results 

As  might  be  expected,  the  modal  patterns  at  low  frequencies  do  not  show 
a  close  relationship  to  mode  shapes  based  upon  classical  boundary  constraints 
and  plate  modal  theory.  For  example,  (b)  of  Figure  10  exhibits  a  modal  pattern 
reminiscent  to  the  first  mode  of  a  fully  clamped  plate  but  (a)  of  the  same  figure 
shows  a  modal  pattern  similar  to  the  fourth  mode  of  a  cantilever  plate.  Obser¬ 
vation  seems  to  confirm  this  latter,  capricious  behavior.  During  this  mode  (79 
CPS)  the  plug-in  (tongue)  end  of  the  board  tended  to  displace  against  the  spring 
contacts  of  the  mating  connector,  in  a  manner  resembling  a  cantilevered 


response, 


At  335  CPS  the  modal  pattern  appears  to  develop  into  a  more  regular  form, 

say  characteristic  of  the  third  mode  of  a  hinged-hinged  plate.  We  explore  this 

possibility  by  examining  the  modal  equation  for  plates.  The  number  of  modes 

of  a  simply  supported  plate  whose  orthogonal  frequencies  (c o  )  are  less  than 

mn 

a  specified  frequency,  cu  ,  is  given  by  (Reference  3): 


N 


(10) 


This  expression  has  a  quadratic  solution  shown  in  Appendix  III. 

Setting,  N  =  3  and  solving  for  the  frequency  (Appendix  III)  we  obtain 
f  =  367  CPS.  This  value  is  in  rough  agreement  with  the  observed  frequency. 
However,  the  results  also  state  that  below  367  CPS,  board  A  contains  three 
mode  patterns.  Examination  of  Figure  10  shows  that  there  are  at  least  five 
modal  patterns  present.  Continuing,  we  set  N  =  1.  Here,  the  frequency  solu¬ 
tion  for  only  one  mode  would  be  an  estimate  of  the  first,  or  fundamental  fre¬ 
quency.  Proceeding  as  before,  we  obtain  f  =  214  CPS.  There  is  no  agreement 
with  the  observed  modal  patterns  —  and,  as  before,  we  note  that  below  214 
CPS  there  are  at  least  three  other  patterns.  Modified  forms  of  Equation  10 
accounting  for  other  boundary  conditions  were  similarly  investigated  —  the 
agreement  was  no  better  and  the  observed  modes  always  greatly  exceeded 
theoretical  predictions. 
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b.  Conclusions 

We  are  now  able  to  account  for  the  peculiar  "bunching"  behavior  of  the 
cumulative  mode  count  at  the  low  frequency  extremity  of  the  curves  of  Figures 
4  and  5.  It  is  evident  from  the  modal  behavior  patterns  of  Figures  10  and  11 
that  the  increasing  modal  density  at  this  portion  of  the  curve  is  due  largely  to 
real  and  variable  boundary  conditions  and  to  the  component  influences  which 
result  in  the  generation  of  more  low  frequency  modes  than  the  predictive  theory 
of  Equation  10  permits. 

It  is  probable  that  most  equipment  configurations  of  this  sort  will  exhibit 
like  behavior  and,  in  the  future,  it  should  prove  profitable  to  subject  other 
plug-in  boards  to  similar  studies. 

Finally,  it  is  instructive  to  note  that  care  should  be  exercised  (when  deter¬ 
mining  modal  densities  via  the  slope  method)  to  insure  that  the  slope  line  is  * 

drawn  above  the  bunched  region  of  the  curve  and  below  the  region  of  modal 

blending;  if  an  accurate  estimate  of  the  modal  density  is  to  be  realized.  V 
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SECTION  IV 

LOSS  AND  COUPLING  FACTORS 


1.  LOSS  FACTOR  MEASURABLES 


Certain  classes  of  loss  factors  must  be  determined  before  the  coupling 
factor  can  be  described.  In  Reference  2  we  note  that  for  a  two  element  energy 
sharing  structure  we  have  two,  interrelated  expressions: 


*7.2  \ 

^2.  +  \ 


^20  ~  % 


V  % 


(ID 


(12) 


Where  rj^  and  77^  are  the  apparent  dissipation  loss  factors  of  the  box  cover 
and  circuit  boards  as  measured  in  their  normal,  interconnected  state,  and  77 
and  Tj2  are  the  true  dissipation  loss  factors,  respectively,  of  the  outer  cover 
as  measured  without  the  presence  of  the  boards,  and  of  the  boards  as  measured 
without  the  presence  of  the  cover. 


From  the  above  expressions,  we  must  obtain  77  so  that  we  can  define 
the  coupling  factor, 


A 


77  +77 

'21  '2 


Selecting  Equation  11,  noting  from  Reference  1  that  77^  N1  =  77^ 
and  solving  for  77^,  we  have: 


_  .  \  \  ~  ^2 

'21  N2 

-  VmT^ 


(13) 


Also,  since  -  n^  (f)  •  Af  and  Ng  -  ng  (f)  •  Af,  we  substitute  in  the  above 
expression  and  obtain: 

^  ^2  ~  *7,0  \ 


'21 


n->  (f ) 


(14) 


v.  -  V.  ~ 


'10  '1  n((f)  ^2 
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We  have  seen  earlier  that  the  total  modal  density  of  the  circuit  boards  is 
very  nearly  equal  to  the  total  modal  density  of  the  cover. 

Allowing  this  assumption  results  in  a  slight  change  of  rj ^  and  an  almost 
imperceptible  influence  on  >  which  is  of  ultimate  importance.  Also,  the 
effect  on  the  prediction  Equation  2  is  small  —  amounting  to  less  than  1.  5  db. 
Thus,  we  set 

n2(f>  , 

n^f)  • 

As  a  consequence,  we  are  permitted  to  set  772a  ~  1?2’  Substituting  these  con¬ 
clusions  into  Equation  14,  we  obtain: 

^2o  ~~  ^ia  ^2a 

77  =  — - (15) 

'21  77-77-77 

'la  'l  '2o 

This  final  expression  now  contains  those  loss  factors  which  can  be  measured 
with  the  exception  of  77^,  which  will  be  discussed  later. 

The  same  recordings  used  for  determining  the  modal  densities  also  pro¬ 
vided  measurements  of  the  Q  and  its  inverse,  77.  The  loss  factors,  determined 
by  the  half  power  point  method,  were  arithmetically  averaged  over  1/3  octave 
intervals.  Measurements  extended  to  approximately  2  KHZ. 

The  loss  factor  for  the  boards  was  measured  for  two  conditions:  a  loss 
factor  obtained  with  the  boards  mounted  in  the  equipment  box  (in  situ)  and  a 
loss  factor  obtained  with  the  board  mounted  in  the  test  fixture  (Figure  12). 

Two  other  plots  were  provided:  one,  (Figure  13)  of  the  in  situ  loss  factor, 
plotted  together  with  an  i72a  of  circuit  boards  measured  by  another  investigator 
(Reference  6)  during  studies  of  an  Atlas  guidance  computer  (the  second  (Fig¬ 
ure  14)  a  plot  of  the  box  cover  mean  apparent  loss  factor,  7?la). 

The  Atlas  computer,  17  ,  exhibits  much  wider  scatter  than  the  equipment 

77  ,  perhaps  because  of  insufficient  averaging  of  modes  within  the  1/ 3  octave 

2ci 

bandwidths.  Nonetheless,  both  loss  factors  are  in  the  same  range  and  seem  to 
follow  much  the  same  inverse  trend  with  frequency.  A  significant  and  useful 
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THIRD  OCTAVE  BAND  CENTER  FREQUENCY  IN  Hz 

Figure  12.  Comparison  of  Circuit  Board  Loss  Factors  When  Mounted  in  Test 
Jig  and  Then  in  Situ 
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Figure  14.  Mean  Apparent  Loss  Factor  ^  of  Equipment  Cover 
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conclusion  arises  from  the  examination  of  Figure  12.  Note  that  the  loss  factor 
of  the  board,  when  mounted  in  the  jig,  seems  to  change  very  little  when  mounted 
in  situ.  It  is  this  observation  that  permits  us  to  set  rj2  ~  17^,  as  described 
earlier. 

2.  OBTAINING  THE  TRUE  LOSS  FACTOR  OF  THE  COVER,  r)1 

Data  had  not  been  obtained  of  the  true  loss  factor  of  the  cover  77.^  however, 
other  investigators  had  recorded  the  reverberation  time  versus  frequency  history 
of  an  empty  box  of  the  same  material,  thickness  and  of  similar  size,  shape,  and 
construction. 

Since  77  =-|p-  where  f  is  the  1/3  octave  center  frequency  and  Tg  is  the 
s 

60  db  decay  time,  we  are  able  to  plot  rj^  and  include  it  with  77^  and  77^. 

This  is  shown  in  Figure  15. 

3.  COMPUTING  T721  AND  fjL^ 

The  three  curves  are  then  represented  as  being  approximately  a  straight 
line,  log-log,  negatively  sloping  function  that  are  asymptotic  on  both  ends  of 
the  frequency  spectrum  (Figure  16).  The  asymptotes  result  from  an  examina¬ 
tion  of  the  trend  of  the  data  and  of  similar  data  from  other  investigators.  Since 
the  measurements  extend  to  only  2000  CPS,  dashed  lines  are  provided  to  indi¬ 
cate  extrapolation  beyond  this  limit  for  values  of  77^  and  17 1&. 

Returning  to  Equation  15,  77^  is  calculated  and  plotted  and  the  results 
are  shown  in  Figure  17.  A  plot  of  an  17^  of  the  Atlas  computer  (Reference  6) 
is  included  for  comparison. 

The  two  coupling  loss  factors  do  not  agree  in  magnitude  and  trend  —  they 
diverge,  noticeably,  above  2  KHZ. 

4.  ENERGY  COUPLING  FACTOR  ^l2 

9 

With  7721  and  77^  defined,  we  plot  ^  as  shown  in  Figure  18.  The  curve 
exhibits  interesting  properties.  Recalling  that  T^  =  T  2  where  T  ^  and  T2 

are  the  modal  vibration  energies  of  the  cover  and  boards  respectively,  we  note 
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2 

that  /i^  has  a  gradual  negative  slope  of  approximately  1.  2  db/octave.  The 
curve  indicates  that  the  printed  circuit  boards  are  behaving  essentially  as  a 
broad-band  vibration  receiver  (that  is,  with  respect  to  their  energy  source, 
the  cover). 
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SECTION  V 

VIBRATION  PREDICTION 


1.  PREDICTION  AND  MEASUREMENT 


Restating  Equation  2  we  have: 

<  S(ou  )2  >  = 


<  S(cu),  > 


n2  (cu)m| 
n,  (cu)m2 


m  is  obtained  by  weighing  the  removed  printed  circuit  boards,  and  amounts 
to  7.  88  pounds.  The  remainder,  m^  is  9.50  pounds.  The  ratio  m^/m^  is  equal 
to  1.  22,  or  1.  7  db  (20  log^.  Since 

n2(tu) 

n,  (cu)  =  1 


Equation  2  becomes: 

<S(^)2>  =  <  Situ),  >  1.22 /i.,2  06) 


We  will  record  our  noise  spectral  values  of  grms  (g)  over  1/3  octave  band- 
widths,  thus  Equation  16  is  expressed  as: 


<  g-(cu)2> 


<g  (tu),  > 


(17) 


or,  in  db  form: 

20  log|0  <  g  (cu  )2>  =  20  >oglo  <g{cu)(>  +  0.9  db  +  10  log  fL,2 

Two  small  high  frequency  shakers  were  attached  to  opposite  sides  of  the  box 
cover  (Figures  19  and  20).  A  random  noise  signal  was  passed  through  aB  &  K 
1612  S/2  spectrum  shaper,  and  from  there  to  the  shaker  amplifiers.  Therefore, 
the  shakers  were  operating  in  parallel  with  respect  to  spectral  shaping. 


Two  accelerometers  (Clevite  2E1)  were  located  on  the  two  opposite  and 
largest  panels  of  the  box.  These  were  designated  panels  one  and  two.  One 
Endevco  2222  accelerometer,  mounted  on  each  of  three  boards  in  the  first 
column  of  boards  (Type  A),  provided  the  response  measurements  for  <g(<u)2>  . 
The  two  cover  accelerometers  provided  a  control  for  spectral  shaping  of  the 
box  cover.  A  curve  was  obtained  to  approximate  a  vibration  response  field 
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Figure  19.  Shaker  Arrangement  and  Accelerometer  Locations  for  Measuring 
Vibration  Level  of  Box  Cover  and  Boards 
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Figure  20.  Shaker  Arrangements  and  Accelerometer  Locations  for  Measuring 
Vibration  Levels  of  Box  Cover 
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resulting  from  coupling  to  an  acoustic  source.  One  accelerometer  was  left  in 
place  after  spectral  shaping  to  provide  a  spectral  control  check,  and  the  re¬ 
maining  accelerometer  provided  the  roving  pickup  points  for  the  box  panel. 
Three  points  per  panel  were  recorded  in  the  1/3  octave  frequency  intervals 
fro  40  Hz  to  16  KHz,  for  a  total  of  18  samples. 

During  each  data  run  on  the  box  cover,  one  location  of  each  of  the  three 
board  accelerometers  was  similarly  recorded.  When  the  measurements  on  the 
three  board  locations  were  completed,  the  accelerometers  were  relocated  to 
the  three  other  boards  located  in  the  adjacent  column  of  boards  (Type  B).  Thus, 
there  were  also  a  total  of  18  readings. 

The  results  of  each  board  were  arithmetically  averaged;  then  the  three 
boards  were  in  turn  averaged.  The  Type  B  column  of  boards  were  similarly 
treated,  and,  finally,  both  board  columns  were  averaged.  This  provided  the 
spatial  response  average,  <g(tu)0>. 

La 

The  same  averaging  procedure  was  applied  to  the  accelerometer  data  from 
the  box  cover  to  obtain  <  gfcu)^  ,  (Figure  21).  As  a  final  step,  the  cover  vi¬ 
bration  field,  <g(cu)1>  ,  was  substituted  into  Equation  17  and  the  predicted 
response  field,  <  g(u))^  >  ,  was  obtained.  This  field  was  compared  to  the  mea¬ 
sured  value.  Both  are  shown  in  Figure  22. 

a.  Results 

The  agreement  between  theory  and  measurement  is  surprisingly  good, 
especially  above  200  Hz.  That  the  spread  between  observed  and  predicted 
values  increases  below  this  region  is  an  expected  result.  This  is  because 
power  flow  between  the  ensembles  of  mode  sets  (in  this  case  represented  by 
1/3  octave  intervals)  requires  that  a  mode  or  modes  exist,  and  rather  uniformly 
so,  in  the  ensemble  interval.  These  conditions  become  less  probable  as  the 
ensemble  bandwidths  decrease  which  in  this  case  with  the  1/3  octave  resolution. 
This  is  in  violation  of  the  funadmental  assumptions  that  permit  use  of  the  power 
flow  equations. 

Starting  approximately  one  octave  above  the  first  modes  of  the  board  and 
box  cover  (approximately  200  Hz)  the  measured  and  predicted  response  values 
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Figure  22.  Comparison  Between  Circuit  Board  Response  Measurements 
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of  the  boards  agree  within  about  3  db  (average)  and  agreement  extends  over 
six  octaves  —  an  impressive  demonstration  of  the  utility  potential  and  extensive 
coverage  capabilities  of  the  energy  method. 

b.  Experimental  ft^2 

If  we  insert  the  measured  values  of  <g(tu)1>  and  <g(cj)2>  into  Equation 
17  we  may  then  solve  for  /x  2,  term  it  experimental,  and  compare  this  coupling 
factor  to  the  previously  derived  form  based  on  a  two  nodal  model.  This  should 
give  us  some  measure  of  the  model  validity,  and  should  provide  us  with  a  use¬ 
ful  comparison  between  experiment  and  theory. 

c.  Discussion 

The  two  coupling  factors  are  shown  in  Figure  23.  Both  curves  agree  in 

o 

trend  but  the  experimental  appears  to  be  higher  everywhere.  Between  500 
Hz  and  5  KHz  both  appear  to  be  on  the  average  within  3  db  of  each  other.  There 
is  a  suggestion  (below  8  KHz)  that  /x^  should  be  somewhat  steeper  in  slope 
than  was  obtained  in  the  derived  form.  This,  in  turn,  suggests  that  should 
also  decrease  somewhat  more  sharply  with  frequency  (Figure  18).  At  this  point, 
it  should  be  noted  that  the  coupling  loss  factor,  Vgi*  a  cumbersome  expres¬ 
sion,  and  is  subject  to  considerable  potential  error. 

Why  the  experimental  coupling  factor  remains  greater  than  /x  derived,  is 
not  clear.  Errors  sttributable  to  the  structural  mass  measurements  and  the 
mode  density  assumptions  do  not  seem  sufficiently  large  to  account  for  a  3  db 
difference.  Thus,  our  suspicion  is  reinforced  that  the  basic  form  of  is  at 
fault.  However,  much  more  detailed  information  concerning  the  nature  and 
extent  of  errors  accumulated  during  the  derivation  of  is  required  before 
this  can  be  firmly  concluded. 

The  most  striking  difference  between  the  two  curves  appears  at  approxi- 
mately  8  KHz.  Here,  /x  measured,  rises  rather  abruptly  and  continues  to 
do  so  up  to  16  KHz.  In  previous  works  (Reference  2)  it  has  been  noted  that  at 
high  frequencies,  vibration  of  the  boards  in  shear  may  result  in  a  higher  mode 
density  than  obtained  for  the  board  during  bending  at  the  lower  frequencies.  It 
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2 

might  be  supposed  that  this  would  account  for  the  behavior  of  since  by 
inspection  of  Equation  2  an  increase  in  would  produce  a  proportionate 

increase  in  the  mean  square  of  the  board  response.  However,  computational 
checks  indicate  that  noticeable  shear  contributions  would  occur  above  approxi¬ 
mately  20  KHz.  Thus,  this  explanation  does  not  appear  satisfactory.  Apart 
from  possible  defects  in  the  assumed  model,  it  may  be  that  torsion  modal 
contributions  from  the  bars  and  stringers  (ignored  in  this  study  in  order  to 
simplify  the  model)  may  account  for  the  observed  increase.  Again,  future 
studies  in  this  area  should  provide  a  more  acceptable  explanation. 

If,  in  review,  we  consider  the  extensiveness  of  the  input  and  response 

2 

measurements  and  weigh  them  against  the  errors  in  the  derivation  of  /i-1  then 

2  x 

it  seems  prudent  to  consider  the  experimental  value  of  as  representing  the 
best  estimate  of  the  coupling  factor  —  at  least  below  8  KHz. 

2 

Accordingly,  a  smoothed  version  of  (Figure  24)  measured,  is  pre¬ 
sented.  The  questionable  region  above  8  KHz  is  indicated  by  a  dashed  line. 

It  will  be  interesting  to  see,  in  the  course  of  future  studies  and  prediction 
efforts  concerning  similar  equipment  structures,  whether  or  not  this  coupling 
factor  will  be  sufficiently  tractable  for  general  use.  In  the  interim,  it  is  so 
proposed. 
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SECTION  VI 
SUMMARY 


1.  CONCLUSIONS 

The  results  of  this  investigation  demonstrate  that  statistical  energy  methods 
may  be  effectively  applied  to  at  least  certain  classes  of  equipment,  consisting 
of  planar  arrays  enclosed  by  a  cover,  and  can  provide  a  useful  prediction  of 
the  internal  vibration  levels  over  an  extensive  frequency  range.  It  is  further 
shown  that  such  complex  equipments  may  be  represented  as  a  two  nodal  sys¬ 
tem,  thereby  simplifying  the  prediction  process.  Other  results  show  that  modal 
densities  of  printed  boards  may  be  determined  from  thin  plate  theory  and  that 
acceptable  engineering  estimates  may  be  realized  by  applying  compensating 
procedures  to  account  for  mass  and  stiffness  effects  of  heavily  laden  boards. 

Also  in  contrast  to  modal  density  estimates,  it  is  demonstrated  that  the  coupling 
loss  factor  is  difficult  to  derive  and  can  normally  be  determined  only  by  indi¬ 
rect  means. 

A  generalized  coupling  factor  is  presented  for  application  to  those 

equipments  having  structural  arrays  and  components  that  are  considered  simi¬ 
lar  to  the  equipment  package  studied  in  this  work.  Finally,  it  is  proper  to  con¬ 
clude  that  as  future  studies  and  measurements  provide  additional  loss  and 
coupling  factors  for  more  diverse  classes  of  equipments;  then  it  will  be  possible 
to  effectively  extend  the  application  of  energy  method  to  more  complex  equip¬ 
ment  and  structures. 

2.  RECOMMENDATIONS 

A  follow-on  study  should  be  conducted  to  determine  the  coupling  factor 
when  the  equipment  box  is  directly  excited  by  an  acoustic  source;  for  example, 
a  reverberant  chamber.  If  this  step  is  successful,  then  it  will  be  possible  to 
predict  the  internal  vibration  of  certain  classes  of  equipment  when  exposed  to 
reverberant  acoustic  sources. 

If  this  phase  is  successful  then,  the  way  may  open  to  simulate  or  reproduce 
the  internal  vibration  levels  of  an  equipment  package  by  coupling  small  high 
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frequency  shakers  to  the  box  surfaces ;  thereby  providing  low  cost  and  flexible 
simulation  of  acoustic  fields.  If  this  technique  proves  feasible  then  a  specifica¬ 
tion  test  method  for  equipments  should  be  synthesized  and  included  in  MIL- 
STD-810(B). 
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APPENDIX  I 

DETERMINATION  OF  MODE  DENSITY  FOR  BOARDS  WITH  COMPONENTS 
AND  STIFFENERS  REMOVED 


Applying  center  loading  to  the  stripped  board,  the  stiffness  k  and  deflection  A  are 
recorded. 


Thus: 


— —  Where  :  — — 

A  A 


k 


Fo  r  board  A 


E  .  (  7.9)3  (10)  (13.1) 
4  (6) 


I  .  14  x  I06  #/in.2 


n(f)  =  no.  of  modes/CPS 
k  =  13.1  #/i  n. 


The  modal  density  is, 
n  (  f  )  =  - -p- 

CL» 

and  CL  =  -n/?"  =  y/'  '  14  *  =  77-5  X  '°3  in. /sec 

I  .  95  1 0 


thus, 

n(f)  =  - — — - =.0182  modes/CPS 

77.5  (  I03)  (  .0625) 


or  , 


I  =  7.9  in. 

b  =  width  =  6  in. 

t  =  thickness  =  1/16  in. 

w  =  weight  of  board  =0.24# 

w 

p  - - 

volume  x  386  in^sec 

.  _ 0.24 _ 

8.5(6)  (  .0625)  (386) 

=  1.95  x  10  4#-sec2/in.4 


n  ( f )  =  18.2  modes/ 1000  CPS 
For  board  B 


(7.  I)3  (  1 6 )3 

4  (7) 

I .5  x  IQ6 
2.22  x  10" 


(29.5)  =  ,  , 


=  82  x  10' 


x  I06  #/in.2 


in. /sec 


p 

Ap  =  board  area  =  51  in. 

C^  =  longitudinal  wave 
velocity  (in. /sec) 

k  =  24.5  #/in. 

1=7.1  in. 

Ap  =  59.5  in.2 
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The  model  density  is, 


n(f)  =  - 59,5  (1,73) -  _  020  modes/CPS 

82  (10*)  (.0625) 


b  =  7  in. 

t  =  1/16  in.  =  .0625  in. 

W  =  .32  # 

p  =  - -2S - 

3.72  (.38  6) 

=  2.22  x  I0-4  #  -  sec2/ in.4 


or 


n (f )  =  20  modes/1000  CPS 


Summarizing: 

The  mode  density  for  board  A  is: 

n  (f)  =  18.  2  modes/1000  CPS 
The  mode  density  for  board  B  is: 

n  (f)  =  20  modes/1000  CPS 
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APPENDIX  II 

MODAL  DENSITY  CORRECTIONS  FOR  INCREASED  MASS  AND  STIFFNESS  OF 
COMPLETE  BOARDS  MOUNTED  IN  SITU 


The  board  is  mounted  in  the  test  jig  and  the  stiffness  k  is  determined  as  in  Appendix  I. 

F 


Dividing: 


Thus: 


Solving  for 


Where: 


k 

k' 


w 


w 

CL 

CL 
E  andE’ 

P  and  p' 


n(f) 

n(f)' 


CT  for  board  A 


=  stiffness  of  board  A  (stripped)  (13. 1  #/in) 

=  stiffness  of  complete  circuit  board  mounted  in  jig  (76  #/in) 

=  weight  of  board  stripped  of  stiffners  and  components  (0.  24#) 

=  weight  of  board  with  components  and  stiffeners  (.  62#) 

=  longitudinal  wave  velocity  of  stripped  board  (77.  5  x  10  in/sec) 
=  pseudo,  or  adjusted  C 

Li 

=  Youngs  Modulus  for  stripped  and  original  boards,  respectively 
=  density  of  stripped  and  original  boards,  respectively 
—  mode  density  of  stripped  board 
=  mode  density  of  complete  (and  mounted)  board 
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For  Board  A 

c'l  *  775  «  '02  - 0#  =  72  2  *  '°’v^ 

cj_  =  77.5  x  I03  (1.5)  =  116  x  I03  in. /sec 


Dividing: 


We  obtain: 


n  (f )  = 


t  C. 


by  n ( f ) '  = 


yr 


t  c, 


n(f )  Cj_ 
n(f )'  CL 

c: 

n(f ) '  = - —  n  (f ) 

C, 


For  Board  B 


n  ( f )’  =  77~3_  (18.2)  =  12.1  modes/1000  CPS 

116 


Where: 


C‘L  =  104 


Cj_  =  104 
Cj_  =  104 


I03  a/  ll6—  Q'32—  *  104  x  I03  v/l.27 

V  24.5  1.19  v 

I03  (1.13)  =  118  x  I03  in. /sec 


(18) 


k  =  24.5  #/in. 
w  =  0.32  # 
k'  =  116  #/ in. 
y»'  =  1.13# 

CL=  104  x  I03  in/sec 
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Substituting  in  Equation  18 

"<f)  =  “l^8~  (20)  =  l7’ 

Summarizing: 


For  board 

A, 

n(f )'  = 

12.  1 

For  board 

B, 

n  ( f  )'  = 

17.7 

modes/1000  CPS 

modes/1000  CPS 
modes  /IOOO  CPS 
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APPENDIX  HI 

DETERMINATION  OF  FREQUENCIES  BELOW  WHICH  A  GIVEN  NUMBER 

OF  MODES  ARE  CONTAINED 


For  a  simply  supported  plate  we  have: 

mn<UJl  = 

Where: 

N  =  number  of  modes 

oj  =  frequency  (rads/sec) 

2 

A  =  area  of  plate  (in  ) 

P  =  perimeter  of  plate  (in) 

c  =  -JojKC  =  wave  velocity  (in/sec) 

K  =  radius  of  gyration  of  plate  section  =  y-yry 

t  =  plate  thickness 

CT  =  longitudinal  wave  velocity 
JLi 

We  set  N  =  3,  and  rearranging  the  above  equation,  we  have: 


Ap  W 
(c)2 


-  P 


(cu) 

(c) 


+  7T 


Where: 


UoT 

( c)2 


-P 


(Cl>) 

(c) 


Ax  +  B„  +  C  = 


-  1 1  7T  =  0 

0 


which  has  the  quadratic  form: 


A  =  A 
B  =-P 
C  =  -  II  7T 

x2  =  — ^ —  *  where  c  =  ^/cuKCL 


Thus 


ki 

2tt 


2 
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Solving  for  the  root 


Where: 


x 


f 

f 


(-b)  ±  y  b2 

-  4AC 

2A 

 p  ±  y(p)2 

-(44  tt)  ( 51) 

2A 

29  ±  y  841 

+  7050 

102 

29  ±  7891 

29  ±  88.8 

102 

102 

KC’  (x2) 

2  V 

333  ( x2)  =  333 

=  367  CPS 

P  =  29 

i  n. 

A  =  51 

.  2 

in. 

C  =  -  I  I  7T 

4AC  =  -7050  in.2 


1.05 

( I.05)2 


yjr  leyir 


c'L=  116  x  I03  in. /sec 


For  the  case  N  =  1,  the  modal  equation  is: 


A 


P 


(u ,) 

-  -  3  v 

(c) 


0 


Proceeding  as  before,  f  =  214  CPS. 
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